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Abstract: Mangrove forests, characterized by high photosynthetic rates and low light 
compensation point, exhibit high Gross Primary Production (GPP), an important component of 
"blue carbon". Accurate estimation of regional GPP and quantification of its limiting factors are 
greatly significant for China to achieve its carbon neutrality goals. In this paper, we estimated the 
GPP of mangrove ecosystems in Hainan Island from 2016 to 2020 based on the Mangrove 
Vegetation Photosynthesis Light Use Efficiency (MVP-LUE) model using Sentinel-2 imagery and 
environmental data; and we also explored the spatiotemporal distribution pattern of GPP and its 
driving mechanism. The results were as follows: (1) During the study period, the annual GPP of 
mangrove ecosystems on Hainan Island showed an increasing trend, with higher GPP in the 
eastern than western regions and northern than southern regions. The mangroves distributed over a 
large area in northeastern Hainan Island dominate the temporal variation patterns at the whole 
island scale. However, distinct differences exist in the temporal dynamics across different regions 
of the mangrove ecosystem in Hainan Island. (2) In terms of the formation mechanism, the spatial 
and temporal distribution pattern of GPP of mangrove ecosystems on Hainan Island was driven by 
a combination of multiple meteorological factors. Seasonally within each year, during the 
transition from the dry season to rainy season, GPP was higher due to lower limitation from air 
temperature (Tair), and the promoting effects of photosynthetically active radiation (PAR) and sea 
surface temperature; in the middle of the dry season, low Tair imposed a serious limitation on GPP, 
but this Tair limitation was weakened with decreasing latitude. In the rainy season, higher cloud 
cover resulted in PAR becoming a limiting factor for GPP. At the end of the paper, we discussed 
the uncertainties and limitations of MVP-LUE, and the spatiotemporal distribution patterns of 
mangrove ecosystem GPP not only be constrained by environmental factors, but species 
composition and forest age structure can also be driving factors. The results of this study provide 
basic data to assess the contribution of regional mangrove forests to the global carbon cycle, and 
theoretical support to reveal the key environmental factors affecting mangrove ecosystem carbon 
dynamics. 

Key words: MVP-LUE model, Sentinel-2, spatiotemporal distribution, limiting factors, mangrove 
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红 树 林 作 为 最 富 碳 的 生态 系统 之 一 , 年 碳 吸 收 速率 为 成 熟 热带 森林 的 2-4 fi, 其 地 下 碳 
储量 可 达 热 带 雨林 的 4-18 fi (Baba etal., 2004) ， 是 全 球 重要 的 碳 汇 (Alongi, 2012) 。 准 
确 估 算 红 树林 总 初级 生产 力 (gross primary production, GPP) 对 于 了 解 红 树林 碳 循环 过 程 及 
评估 其 对 全 球 碳 循环 的 贡献 非常 重要 ; 此 外 ， 深 入 探究 红 树 林 生 态 系 统 GPP 对 环境 因素 的 
响应 不 仅 有 助 于 增强 模型 预测 的 精度 , 而 且 有 利于 掌握 气候 变化 对 于 红 树 林 生 态 系 统 生产 力 
的 影响 。 

然而 ， 红 树林 分 布 于 海岸 潮 间 带 ， 受 到 陆 海 协 同 作用 影响 ， 估 算 其 GPP 不 仅 需 要 考虑 
典型 陆地 生态 系统 的 环境 影响 因子 [如 空气 温度 Cair temperature, Tair) ~ WAKAAE (vapor 
pressure deficit, VPD) 、 光 合 有 效 辐射 (photosynthetically active radiation; PAR) 等 (Sun et 
al., 2019) ]， 还 应 考虑 到 海洋 水 面 温度 (sea surface temperature, SST) 和 海水 盐 度 (sea water 
salinity, SAL) 等 生境 特性 影响 因子 (Wu et al., 2023) 。 学 者 已 经 取得 了 丰富 的 关于 环境 因 
子 对 红 树 林 GPP 影响 的 科学 发 现 (Rodda et al., 2022), 最 优 光 合作 用 Tai 为 28~32 "CCGilman 
et al., 2008) ， 过 高 或 过 低 的 Tar 均 会 降低 光合 酶 促 反 应 速率 ， 同 时 ， 高 温 也 会 降低 气孔 导 


度 以 应 对 高 蒸 散 发 需求 ， 进 而 限制 GPP (Niu et al., 2012) . VPD 也 是 蒸 散 发 需求 的 表征 因 
子 ， 较 高 的 VPD 将 导致 红 树林 气孔 闭合 (Barr, 2005; Alongi, 2009) ， 限 制 光 合作 用 。 红 本 
林 碳 的 同化 在 较 低 水 平 上 受到 太阳 辐射 的 显著 影响 ， 而 高 PAR 会 给 冠 层 带 来 过 多 的 光 吸 收 
和 热量 , 降低 冠 层 的 光 能 利用 效率 (light use efficiency, LUE ) 导 致 GPP 降低 (Cui et al., 2018)。 
此 外 , 受到 周期 性 潮汐 作用 , 尽管 红 树林 已 经 发 展 出 特殊 的 结构 或 组 织 来 适应 这 种 苛刻 的 环 
境 (Cui et al., 2018) , 1E SST 和 SAL 这 两 个 环境 因子 仍然 是 影响 红 树 林 GPP 的 关键 (Barr 
et al., 2013; Zheng & Takeuchi, 2022) . SAL 在 最 佳 范围 内 (5~15 pp) 对 红 树 生长 和 光合 生 
产 上 具有 强烈 的 正 向 影响 (Liu & Lai, 2019) ，SAL 过 高 则 导致 了 红 树 林 的 根系 环境 中 的 负 渗 
透 压 (Barr et al., 2010; Noor et al., 2015) ， 限 制 水 分 供应 ， 抑 制 光 合作 用 。 潮 汐 淹 没 也 会 改 
变 土壤 温度 , 对 红 树 植物 根部 的 生理 代谢 和 叶 部 的 光合 作用 产生 影响 , 土 温 过 低 会 冻结 红 树 
林 的 木质 部 ， 抑 制 水 分 输送 ， 限 制 光合 活动 (Zheng & Takeuchi, 2022) ; 土 温 过 高 根部 呼吸 
作用 增强 ,大 量 脱落 酸 经 木质 部 输送 至 地 上 部 ， 降 低 气 孔 导 度 ， 抑 制 光合 速率 。 潮 汐 作 用 使 
得 红 树 林地 区 土 温 难 以 测量 , 但 土 温 与 SST 相关 性 良好 , 以 SST 表征 土 温 (Quisthoudt et al., 
2012; Barr et al., 2013) 便 给 红 树 的 相关 研究 带 来 了 便利 。 同 时 ，SST 决定 着 红 树 分 布 的 纬度 
EPR CQuisthoudt et al., 2012; Ximenes et al., 2018) ， 并 且 多 项 研究 也 表明 SST 是 红 树 植物 
生长 的 关键 气候 调节 因子 (Brienen et al., 2010; Ximenes et al., 2018) ， 在 准确 估算 红 树 林 
GPP 方面 表现 出 巨大 的 潜力 。 
在 GPP 的 估算 手段 方面 ， 传 统 的 观测 方法 中 静态 箱 法 灵敏 度 高 、 价 格 便宜 ， 但 箱 体 本 
身 会 扰动 观测 环境 且 难 以 获得 连续 的 时 间 和 空间 数据 ( 宋 霞 等 ， 2003) ， 人 力 成 本 高 。 涡 度 
相关 法 (eddy covariance, EC) 可 以 获取 长 期 连续 的 观测 数据 ， 是 估算 站 点 尺度 的 红 树 林 生 
态 系统 GPP 的 最 优 方式 ， 为 区 域 尺度 GPP 估算 提供 了 校准 参考 依据 ， 可 依旧 面临 通 量 点 位 
少 、 分 布 不 均 且 空间 不 连续 等 劣势 (Liu & Lai, 2019) 。 传 统 估 算 模 型 中 统计 模型 估算 结果 
误差 大 (Yoshikai et aL, 2022) ; 过 程 模型 模型 结构 过 于 复杂 ， 模 型 间 的 差异 甚至 高 于 估算 
的 不 确定 性 本 身 (Sun etal., 2019) ， 并 且 基 于 典型 陆地 生态 系统 所 获得 的 过 程 参数 是 否 适 
用 于 红 树 林 生 态 系统 还 有 待 验证 CYoshikai etal., 2022) 。 随 着 遥感 (remote sensing, RS) 
技术 的 发 展 , 基于 LUE 理论 的 模型 逐渐 成 为 了 区 域 尺度 GPP 估算 的 主要 方法 (Song, 2013 )。 
LUE 模型 由 Monteith(1972 ) 首 次 提出 , 现 已 开发 了 一 系列 广泛 应 用 的 模型 (Sun et al., 2019), 
对 于 评估 生态 系统 GPP 具有 重要 的 实践 意义 。 但 受 限于 红 树 林 独 特 的 生境 ， 且 红 树 通 量 塔 
测量 数量 不 足 CPastorello et al., 2020) ， 准 确 估算 红 树 林 生 态 系统 的 GPP 具有 一 定 的 难度 。 
随 着 Barr 等 首先 根据 包含 海洋 水 文 的 环境 因子 对 红 树 林 LUE. 的 影响 , 建立 了 针对 红 树 林 生 
态 系 统 的 LUE 模型 ， 并 估算 了 在 南 佛罗里达 大 沼泽 地 国家 公园 红 树 林 生 态 系 统 GPP， 且 表 
现 良 好 (Barr et al., 2013) ; Zheng 等 人 在 Barr 模型 的 基础 上 量化 了 SST 和 SAL 对 红 树 林 
生态 系统 GPP 的 影响 ， 改 进 了 红 树 林 植 被 光合 作用 - 光 能 利用 效率 Cmangrove vegetation 
photosynthesis light use efficiency; MVP-LUE) 模型 ， 并 估算 了 中 国 沿 岸 18" 一 28" N 范围 的 
红 树 林 生 态 系统 GPP， 其 准确 度 远 高 于 MODIS 等 遥感 产品 (Zheng & Takeuchi, 2022) ; 
Wu 等 (2023 ) 则 利用 改进 的 MVP-LUE 就 台风 影响 后 的 红 树 林 生 态 系统 GPP 恢复 情况 开展 
分 析 ， 也 得 到 了 MVP-LUE 精度 良好 的 结论 ; 这 都 为 开展 区 域 尺 度 红 树林 生态 系统 GPP 的 
估算 奠定 了 理论 与 技术 基础 。 

准确 估算 红 树 林 生 态 系统 GPP 并 解释 其 时 空 分 布 格 局 对 于 了 解 红 树 林 生 态 系 统 碳 过 程 
具有 重要 意义 。 海 南 岛 红 树林 湿地 资源 丰富 ， 红 树 植物 资源 占 全 国 总 数 的 95%， 也 是 我 国 
红 树 林 分 布 最 集中 、 保 存 最 完整 的 地 区 〈 王 鸿 平 ，2018) ， 是 开展 红 树 林 GPP 时 空 分 布 格 
局 研究 的 理想 场地 。 本 研究 以 海南 岛 红 树林 生态 系统 作为 研究 对 象 ， 利 用 哨兵 -2 影像 与 再 
分 析 气 象 数 据 ， 基 于 海南 岛 本 地 化 的 MVP-LUE 模型 ， 估 算 2016—2020 年 海南 岛 红 树林 生 
态 系统 GPP， 探 讨 海 南 岛 红 树林 生态 系统 GPP 时 空 分 布 格局 及 其 形成 机 制 。 拟 通过 对 上 述 
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号 ”高 分 辩 率 (0.5m 2 
调查 验证 ， 得 到 海南 岛 红 树林 


1.2.2 卫星 数据 


本 而 


究 中 使 用 了 2015 年 1 


BE ZI 
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2310 m, HERI ESEH 
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月 至 2022 年 12 


乡 像 ， 


星 图 


像 ， 


的 Sentinel-2 影 
并 通过 哨兵 数据 应 用 平台 计生 
晶 射 吸收 系数 (fraction of absorbed photosynthetically active radiation，fAPAR) ， 空 间 分 
合成 法 (Holben, 1986) 获得 逐 月 fAPAR. 
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Distribution map of mangrove forests in Hainan Island 


树林 分 布 的 大 致 范围 ， 利 用 2022 年 的 “吉林 一 号 ” 
空间 分 辨 率 ) 多 光谱 卫 
分 布 现 状 数据 。 


通过 目 视 解 译 结 合 实地 


像 ， 包 括 Sentinel-2A 和 
吸收 光合 有 效 


1.2.3 环境 数据 
本 研究 使 用 的 环境 数据 主要 包括 气象 数据 和 海洋 水 文 数据 , 包括 Tar VPD, PAR, SST 
和 SAL. 为 与 fAPAR 的 空间 分 辨 率 一 致 ， 采 用 双 线 性 插值 法 将 环境 因子 数据 重 采 样 至 10 m 
分 辩 率 。 环 境 因子 数据 集 详细 参数 见 表 1. 
表 1 环境 因子 数据 及 来 源 


Table 1 Environmental factor data and sources 


LEES 
数据 时 间 Resolution 产品 名 称 数据 源 
Data Time 空间 时 间 Product name Data source 


Spatial ^ Temporal 


ZH PAR 产品 http://www.glass.umd.edu/i 
光合 有 效 辐射 PAR 2016—2020 0.05° 
Monthly PAR product ndex.html 
https://cds.climate.copernic 
: 每 日 ERAS 再 分 析 数 据 
空气 温度 Tair 2016 一 2020 0.25° us.eu/cdsapp#!/search?type 
Daily ERAS reanalysis data 
=dataset 
https://cds.climate.copernic 
. . 每 日 ERAS 再 分 析 数 据 
海水 表面 温度 SST 2016—2020 0.08° , l us.eu/cdsapp#!/search?type 
Daily ERAS reanalysis data 
—dataset 
2 混合 坐标 海洋 模式 , 
每 日 n a https://code.earthengine.go 
海水 表面 盐 度 SAL  2016—2020 0.08? SAL 数据 集 
Daily ogle.com/ 
HYCOM SAL dataset 


XX H TerraClimate 数据 自 


YT 


https://code.earthengine.go 


饱和 水 汽 压 差 VPD 2016—2020 。 ~4 km 


Monthly TerraClimate dataset ogle.com/ 


1.3” 红 树 植 物 光合 作用 - 光 能 利用 效率 模型 (MVP-LUE) 

本 研究 采用 Wu 等 人 改进 的 MVP-LUE 模型 (Wu et al., 2023) ， 该 模型 是 在 Zheng 和 
Takeuchi (2022) 将 Tars SST, PAR, VPD 和 SAL 作为 影响 红 树林 LUE 的 环境 因素 且 定 义 
相应 的 缩放 因子 建立 的 MVP-LUE 模型 基础 上 , 结合 海南 岛 所 处 纬度 范围 , 确定 了 海南 岛 红 
树林 生态 系统 本 地 化 参数 的 MVP-LUE。MVP-LUE 模型 被 用 来 估算 了 中 国 189-28? N 红 树 
林 的 GPP (Zheng & Takeuchi, 2022) , 并 在 估算 修 州 湾 红 树林 GPP 中 取得 了 较 好 的 结果 (Wu 
et al., 2023) 。MVP-LUE 模型 基础 算法 如 下 : 

GPP = PAR x fAPAR x LUE (1) 

式 中 : GPP 表示 总 初级 生产 力 ; PAR 表示 光合 有 效 辐射 ， fJAPAR 表示 光合 有 效 辐射 吸 
收 系数 ， LUE 则 根据 环境 压力 对 最 大 光 能 利用 效率 LUE 的 缩放 情况 进行 量化 。 量 化 公式 
如 下 : 


"i 


LUE = LUE, xT, san XVPD, uu X SST au XSAL, uu XPAR a Q) 


ir scalar scalar scalar 


式 中 : LUEmax 为 最 大 LUE, Tair scalar、VPDscalar、SSTscalar、SALscalar 和 PARscalar 分 别 为 Tair、 
VPD, SST, SAL 和 PAR 对 红 树 林 LUE 影响 的 限制 因子 ， 模 型 各 参数 与 缩放 因子 的 计算 方 
式 详 见 Wu 等 (2023) 。 

1.4 Theil-Sen 与 Mann-Kendall 趋势 分 析 

非 参 数 统计 的 Theil-Sen 趋势 分 析 和 非 参 数 Mann-Kendall 检验 方法 结合 , 能 很 好 地 应 用 
于 长 时 间 序 列 数据 的 趋势 分 析 中 且 目 前 已 得 到 广泛 应 用 (Sen, 1968; Tucker et al., 1991) , 
中 ，Theil-Sen 趋势 分 析 对 测量 误差 、 离 群 数据 的 敏感 性 较 低 ， 可 以 减少 数据 异常 值 的 影响 ; 


而 Mann-Kendall 是 一 种 非 参数 统计 检验 方法 ， 该 方法 不 要 求 测量 数据 遵循 正 态 分 布 ， 具 有 


抵抗 数据 误差 的 能 力 CMilich et al., 2000) 。 
趋势 分 析 方 法 进一步 分 析 海 南 岛 红 树林 生 
1.5 研究 流程 

本 研究 技术 路 线 如 图 2 所 示 : 
岛 红 树林 分 布 ; 
算 逐 月 环境 压 
统 GPP; (452 
海南 岛 红 树林 生 ; 


ÉH 


会 制 海南 岛 红 树林 生态 
态 系统 月 均 GPP 时 空 


CD 利用 高 分 辨 率 卫 星 

(2) 提取 海南 岛 红 树林 分 布 区 的 环境 因子 数据 ， 并 进行 逐 月 合成 ; 
力 的 缩放 因子 ， 利 用 MVP-LUE EX £55 
系统 GPP 及 各 环境 因子 缩放 值 的 时 空 
分 布 格局 形成 机 制 。 


因此 ， 本 研究 使 用 Theil-Sen 结合 Mann-Kendall 
态 系统 GPP 在 2016—2020 年 的 整体 变化 趋势 。 
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GPP = LUE x fAPAR X PAR 
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海南 岛 红 树林 生态 系统 GPP 时 空 分 布 格局 形成 机 制 分 析 
Analysis on the driving mechanisms of spatiotemporal 
pattern of GPP of mangrove ecosystem in Hainan Island 


图 2 技术 路 线 
Fig.2 Technical flowchart 


2 结果 与 分 析 


2.1 海南 岛 红 树林 生态 系统 GPP 时 间 分 布 特征 


2016 至 2020 年 海南 岛 红 树林 生态 系 
ME dS mt 


年 际 和 年 内 GPP 分 布 如 图 3 所 示 。 海 南 岛 红 树 


年 GPP 在 2016 年 至 2019 年 持续 增加 ， 由 (1 305.35 土 649.46) g C-m?-year'! 提高 


至 (1 591.45 土 747.37) g C-m?-year!, 2019 至 2020 年 有 所 降低 ， 但 整体 上 呈现 出 增长 趋势 ， 
多 年 均值 为 (1 448.58 土 679.06) g Cm?'year1。 海 南 岛 整体 的 多 年 月 均 GPP 呈现 出 年 内 单 峰 特 
征 ，“ 旱 一 雨 ” 过 渡 期 较 高 ， 峰 值 出 现在 4 月 [(184.74 土 99.38) g Cm>:month-!]， 而 旱季 中 期 
最 低 ， 最 低 值 出 现在 1 月 [(77.75 土 42.51) g Cm?-month1]， 月 均 GPP 为 (120.72 土 30) g 


C:m*-month-!。 
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图 3 海南 岛 红 树林 生态 系统 年 际 〈A) 和 年 内 (B) GPP 分 布 


Fig.3 Interannual (A) and intra-annual (B) GPP for mangrove ecosystems in Hainan Island 
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2.2 海南 岛 红 树林 生态 系统 GPP 空间 分 布 特征 
2.2. 年 际 GPP 空间 分 布 与 趋势 分 析 

海南 岛 红 树林 生态 系统 年 均 GPP 的 空间 分 布 ( 图 4: AO 整体 呈现 “ 东 高 西 低 ， 北 高 南 
低 ” 的 空间 分 布 格局 ， 高 值 主要 集中 在 东部 及 东北 部 〈 东 寨 港 、 清 澜 港 ) 。 空 间 上 ， 通 过 
Theil-Sen 值 趋势 分 析 和 Mann-Kendall 显著 性 检验 可 以 发 现 , 海南 岛 80.96% 面 积 的 红 树 林 生 
态 系 统 GPP 呈现 增加 趋势 ， eo ei are 
GPP 降低 趋势 ， 主 要 集中 在 岛 西部 的 东方 与 北部 的 新 盈 湾 等 地 区 ， 也 无 法 通过 显著 性 检验 。 
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图 4 海南 岛 红 树林 生态 系统 年 均 GPP 空间 分 布 〈(A) 和 变化 趋势 (B) 图 
Fig.4 Annual average GPP spatial distribution (A) and change trend (B) maps for mangrove 


ecosystems in Hainan Island 


2.2.2. ”海南 岛 红 树林 月 尺度 GPP 时 空 分 布 格局 
海南 岛 红 树林 多 年 月 均 GPP 时 空间 分 布 如 图 5 所 示 。 海 南 岛 北 缘由 东 向 西 年 内 GPP XX 
峰 特 征 逐 渐 明 


X. 分 别 出 现 在 “ 雨 心 旱 ” 两 季 的 过 


寺 渡 期 ; 海南 岛 东 部 及 东北 部 红 树 林 年 内 单 
峰 特 征 显著 ， 峰 值 出 现在 旱季 向 雨季 转换 的 4 月; 而 南部 地 区 ， 整 体 上 GPP 较 低 ， 并 且 呈 


现 出 了 旱季 高 于 雨季 的 特征 。 从 空间 上 看 年 内 GPP 分 配 ， 海 南 岛 各 区 域 红 树林 生态 系统 逐 
月 GPP 变化 特征 差异 明显 ， 岛 西部 、 北 部 和 东部 红 树 林 生 态 系 统 GPP 在 月 均 尺 度 上 呈现 较 
大 的 波动 ; 南部 的 红 树 林 生 态 系统 逐 月 GPP 变化 相对 较为 稳定 。 
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图 5 海南 岛 典 型 红 树林 区 域 逐 月 GPP 时 间 CAO 和 空间 B) 分 布 图 
Fig. 5 Monthly GPP temporal (A) and spatial (B) distribution maps for typical mangrove regions 


in Hainan Island 


2.3. ”时空 分 布 格局 的 形成 机 制 分 析 


本 研究 利用 MVP-LUE 模型 的 环境 限 


VPDscalar〉 表征 红 树 林 生 态 系 统 GPP 的 环境 压力 ， 各 


各 区 域 不 同时 间 段 的 主导 限制 天 


区 SST 主导 雨季 后 期 ， 在 北部 和 西 


较 高 ， 在 东部 地 


强 于 东部 与 南部 。 与 SST ARIS, Ta 同样 作为 环境 温度 压力 ， 对 各 区 域 纪 
分 布 特征 较为 一 致 , 旦 双 峰 ,时 


随 着 纬度 的 降低 ， 红 树林 GPP 受 Ta FER S 


期 (12 


spl 
至 SE 2 


部 SST 则 会 主 


fi] 29 


HAI (PARscalar、SALscalar、 SSTscalar、 Tair scalars 

素 的 年 内 分 布 特征 如 图 6 所 示 ， 
素 不 同 。 全 年 SAL 的 限制 变 幅 较 低 ; 而 SST 限制 的 变 幅 则 
导 整 个 雨季 ， 并 且 限 制 也 
[ 树 GPP 显著 的 年 内 


季 与 雨季 中 期 的 制约 性 最 强 , 雨 旱 两 季 过 渡 期 的 限制 则 较 低 ; 


东部 地 


主导 着 大 部 分 地 


A 


ru 


限制 ， 鸟 北 地 区 时 月 


过 渡 期 、 南 部 早 季 中 期 至 末期 也 主要 受到 PAR 的 制约 。 


判 程 度 减弱 ; 主导 东部 、 东 北部 及 北部 地 区 旱季 中 
左右 ) 红 树林 生态 系统 GPP。 水 分 压力 (VPDscaar) 与 辐射 资源 限制 
(PARscaar) 表现 出 一 致 的 年 内 单 峰 变 化 特征 ， 雨 季 的 限制 强 于 旱季 ;4 
系统 GPP 均 未 受到 水 分 压力 主导 ， 但 辐射 限制 却 
区 清 澜 港 与 琼 海 纪 


岛 各 处 红 树林 生态 
区 红 树 林 生 态 系 统 GPP; 
[树林 生态 系统 由 早 季 向 雨季 过 渡 直 至 雨季 中 期 均 受 到 PAR. 的 强烈 


aiot ^ F 210 全 
gore " is glor etes 35 g "E " ed 
9 * . 2 * E x - v 
总 09| a d M » * " ELL x 1 * * 5 昌 09 Feet » b: 
oD ew. cre gs 3. oD ree ne v. * "T ! -ge z 
[EE i». n 信 三 08 v 全 .ES08 a 
m * 本 号 a gr = * LI 
Z Eor Bs. K Eo -je EST ii 
Ei 8 s 
m gos j m pos mE 
Li Bos HB [3 os ENDE 
g^ í £04 E04 
go 中 
"ui b [:] d ui e 
03 ar 03 03 
—————— — S 上 n á 
1 3 56789101 0234567 9 10 1 12 1 2 3/4 5 6 7 8 9 10 II 12 pr 3 5678 10 11 1 
月 Month 月 Month N / 月 Month 月 Month 
alo] e— e. Kc QM tg å eh 
g Y $ » 
Ruban A Ni uet. Bp Wem" ta 
a t> ese e Š Wes o v9 e. 
E * A x 
E A JT p EU - — 
E Eo7 —— ED Eo: - . 
2E | m 
E so6r E06 
B5 í mE 
区 Eosr 长 Hall 
e e 
S04 E04 
E a E f 
Hos BosL 
——— — —. a " — 
1234567829 10H1Hltl $4 85 TY $ 9 IU 114 
H Month H Month 
10 - "ELI hé a ^ 
" ese, om 124 E a mate i 
$ wf $ d» t 
£o» * Ah Zo S 
9 M E] Y p 
& EE g pooo pigos E ? . 和 m 
0.8 * í o B08 * 
p . j 由 三 * Dr 
a Zo -iee z 3 : 4 
gE- EZ 
EZ ES 6 
LE [T 
i$ 895 & 8 05 
Em 1 ESI 
三 & 
Hi [5] g 
EIE 3 
A -— 
1 34567829101 1 2345677829 101112 
月 Month N 月 Month 
10 . p10 ALO Kii H ^ 
BC etes pt cw MT se |. ya ed 
$ * e 入 多 E x * $N 
909. v .» xv v art v *»— 3 Y A09 Yo BoF a re ax 
o toene pog 9o e. | m (| me e t] toge r1 ——— ——— I-A e m€"—— X € 
y E08 . pE 08 "- os 1 4 
ng On e— y EE | ad mag- dee, S 入 z . gr 
= E07 p-e x A E Eo nid "n Š E07 oeat 
RE EX um 
= a. í rim = 
FE i: ES 06 & S os 
^o 9 2 
& Eos | Eos Bos 
[3 £ e 
gua , EM n EXE 
E k J E os 1 P as h 
123456278 9101 102 1234567 8 9 101112 123454672829 101 12 2345267 8 9 101112 
H Month 月 Month H Month 月 Month 


图 6 限制 因子 对 海南 岛 红 树林 生态 系统 月 均 GPP 时 空 分 布 格局 的 影响 


Fig.6 Influence of limiting factors on spatiotemporal distribution patterns of monthly average 


GPP in mangrove ecosystems of Hainan Island 


3 讨论 


3.1 模型 不 确定 性 


本 研究 利用 Sentinel-2 影像 和 环境 因子 数据 ， 基 于 MVP-LUE 模型 (Zheng & Takeuchi, 


2022) 估算 了 10m 4 EX F 2016—2020 年 海南 岛 红 树林 生态 系统 GPP。 通 过 对 比 现 有 下 


究 发 现 该 MVP-LUE 模型 在 估算 海南 岛 红 树林 生态 系统 GPP 方面 具有 很 高 的 适用 性 。 其中， 


Zheng 等 使 用 该 MVP-LUE 模型 估算 的 中 国 红 树 林 GPP 其 研究 区 域 已 包含 海南 岛 (Zheng & 
Takeuchi, 2022) ， 估 计 值 的 均 方 根 误差 (Root Mean Square Error, RMSE) 为 0.005 1~0.006 


8 g C m?-d^! 与 实测 值 趋势 接近 


; 表现 显著 高 于 MOD17 (RMSE=0.021 7~0.022 0 g C-m?-d^); 


Wu 等 (2023) 使 用 该 MVP-LUE 模型 则 直接 估算 了 海南 岛 任 州 湾 红 树林 GPP， 其 估计 值 与 


MODIS GPP 产品 经 过 交叉 验证 得 到 其 RMSE 分 别 为 (0.61 和 1.62 g C-m?d) ， 也 更 相对 


准确 。 此 外 , 该 模型 相 较 于 现 有 


的 模型 在 模拟 性 能 的 提高 上 具有 更 大 的 优势 , 如 杨 匡 翔 (2020) 


将 SAL 参数 化 耦合 进 光 能 利 上 


日 率 模型 后 的 LUE sa 模型 所 估算 的 云 雷 站 和 雷州 站 的 红 树 林 


GPP (RMSE=0.99 和 0.88 g C-m?-d) 相 比 于 VPM 模型 (RMSE=1.08 和 1.14 g C-m?-d^) 
估算 的 结果 在 模型 性 能 的 提升 方面 较 明 显 。 因 此 ，MVP-LUE 模型 能 够 很 好 地 用 于 估算 海南 
岛 红 树林 生态 系统 GPP ( 表 2) 。 


然而 ， 基 于 MVP-LUE 模型 估算 的 红 树林 生态 系统 GPP 还 存在 一 定 的 不 确定 性 ， 其 不 


确定 性 主要 来 源 于 : (1) MVP-LUE 模型 的 参数 是 针对 红 树 林 生 态 系统 ， 而 不 是 某 一 种 具 


体 的 红 树 林 植 物种 类 。 不 同 红 树林 植物 种 类 与 环境 关系 的 差异 使 红 树 林 具 有 不 同 的 最 适 生长 


环境 条 件 (Call et al., 2019) ， 


比如 老鼠 各 不 耐 盐 〈( 黄 源 欣 等 ，2023) ， 而 同 为 真 红 树 植物 


的 白骨 壤 则 耐 盐 和 耐水 淹 的 能 力 强 , 是 中 国 红 树 林 中 抗 盐 性 最 强 的 红 树林 先锋 树种 ( 杨 锋 等 ， 
2021) 。 不 同 红 树林 植物 种 类 耐 受 的 环境 压力 (生态 幅 ) 差异 会 带 来 一 定 的 不 确定 性 。 其 次 ， 
植物 在 不 同 发 育 阶段 的 LUEaax 也 存在 差异 〈Coops etal., 2010; 高 德 新 等 ，2021) ， 因 此 ， 


时 龄 对 参数 设置 所 带 来 的 不 确 


定性 也 不 可 忽视 ， 较 好 的 一 点 是 MVP-LUE 模型 的 参数 中 含 


fAPAR， 与 归 一 化 植被 指数 (normalized difference vegetation index, NDVI) 具有 和 良好 的 相 
KRAKA (Yuan et al., 2019) ， 能 够 表征 林 龄 、 营 养 、 冠 层 结构 等 综合 性 状 ， 一 定 程度 上 会 降 


低 模型 不 确定 性 。 (2) 不 同 环 
在 一 定 程度 上 提高 了 温度 条 件 


境 变 量 之 间 存 在 共 线 性 作用 。 例 如 Ti 和 SST 限制 相互 关联 ， 
的 约束 (Call et al., 2019) ， 因 此 ， 不 同 环 境 变 量 之 间 的 交互 


性 不 明确 会 增加 模型 估算 的 不 确定 性 。 (3) 遥感 数据 相对 较 低 的 时 间 分 辨 率 〈 颜 秀 花 等 ， 


2019; R, 2020; RIE, 2023) 以 及 再 分 析 数 据 较 低 的 空间 分 辨 率 均 会 影响 估算 的 


准确 性 。 
本 研究 中 未 涉及 降水 影响 
70%~90% (Cui et al., 2022) ， 


， 主 要 考虑 到 海南 岛 雨 季 降 雨量 较 高 ， 可 达 全 年 总 降雨 量 的 
其 中 以 热 雷 雨 和 台风 雨 为 主 ， 二 者 降水 量 可 达 90% 以 上 。 热 


雷雨 及 台风 雨 短 时 、 快 速 、 大 量 ， 大 部 分 形成 径流 无 法 被 生态 系统 所 利用 ， 因 此 在 月 尺度 上 
将 降水 量 作为 环境 压力 会 带 来 较 大 的 不 确定 性 。 但 水 资源 作为 生态 系统 重要 的 环境 因子 , 是 
需要 被 考虑 的 , 本 研究 中 采用 VPD 与 SAL 并 结合 PAR 进行 代理 , 水 资源 充足 时 VPD 较 低 ， 


SAL 可 以 衡量 红 树 林 生 态 系 统 


降低 降水 量 对 于 MVP-LUE 模型 的 不 确定 性 。 
表 2 本 研究 所 采用 MVP-LUE 模型 的 精度 分 析 


的 淡水 输入 量 ，PAR 较 高 则 表征 云 覆 盖 低 降水 量 不 足 。 以 此 


Table 2 Accuracy analysis of the MVP-LUE model adopted in this study 


MODIS 产品 或 其 他 LUE 模型 估 


本 研究 的 均 方 根 误差 、 
研究 区 算 的 均 方 根 误差 来 源 
RMSE in the study 
Study area RMSE of MODIS products or References 
(g Cm?) 


other LUE models 


(g Cm?d!) 


Zheng & 
中 国 红 树林 
0.005 1~0.006 8 0.021 7~0.022 0 Takeuchi, 
Mangrove forests in China 
2022 
ir f e SEE RE 
Wu et al., 
Mangrove forests in Danzhou 0.61 1.62 
2023 
Bay, Hainan Island 
福建 省 云 直 站 红 树 林 , 
. AR, 
Mangrove forests at Yunxiao 0.99 1.08 
2020 
Station, Fujian Province 
广东 省 雷州 站 红 树 林 
| [LEN 
Mangrove forests at Leizhou 0.88 1.14 
2020 
Station, Guangdong Province 
3.2. 时空 分 布 格局 的 形成 机 制 


海南 岛 整体 的 逐 月 GPP 变化 特征 与 东部 、 东 北部 乡 


与 东北 部 红 树 林 分 布 面积 / 


> É 


积 权重 较 高 。 在 旱季 


及 南海 北部 的 海 温 偏 低 
12 至 次 年 2 月 气温 偏 低 〈( 朱 


要 控制 


Ed Ej 
TERS 


且 受 海洋 热力 强迫 作用 , 使 得 冷 


[树林 
期 (1 


ZH GPP 相似 ， 主 要 | 


于 东部 


至 次 年 2 月 ) 我 


20180 ， 低 温 被 


为 适合 红 树 林 生 长 的 温度 条 件 是 : 最 冷 月 均 温 度 高 于 20 


型 温度 ， 低 温 可 能 会 抑制 红 树 的 光合 速率 ， 导 致 了 Tar 是 海南 岛 纪 


S 


pr 


LIUNS: 


pA 


°C, H. 


至 次 年 2 


子 ， 但 由 于 未 完全 进入 雨季 ， 


主要 的 环境 限制 因子 ,但 是 随 着 纬度 的 降低 而 Tai 逐渐 升 
到 Tar 的 限制 有 所 减弱 。 到 了 旱季 末期 和 雨季 初期 (3 一 5 月 ) PAR 逐渐 成 为 主要 限制 因 
云 履 盖 度 尚未 达到 峰值 ， 因 


国 东部 近海 
向 海南 岛 ， 
泛 认 为 是 全 球 红 树 林 纬 度 范围 的 主 
因素 (Lugo & Zucca, 1977; Tomlinson, 1986; Duke et al., 1998) , Walsh 等 (1974) iA 
季节 温差 不 超过 5 °C 的 热 
[树林 生态 系统 GPP 在 12 
高 ， 使 得 海南 咏 南 


使 海南 岛 


t 


F^ 


TH 


部 受 


此 PAR 限制 3 


不 强 ; 同时 温度 的 


回升 ， 逐渐 接近 最 优 光 合 温度 ,提高 了 光合 酶 促 反 应 速率 ; 渐 入 雨季 湿度 同时 也 不 是 限制 红 


树 光 合生 产 的 因子 ; 综 上 ， 海 南 岛 纪 
峰值 。 而 后 
红 树 3 
呈现 | 


季 的 降水 湿度 也 较 大 ， 
和 东部 地 

背风 坡 ， 降 水 量 少 ， 受 
WRK, 雨季 期 间 湿度 


西北 季风 强 弱 


西北 向 东南 递减 空间 分 布 格局 ， 因 而 海南 
递减 。 然 而 ，SAL 在 不 同月 份 下 和 位 置 上 都 较为 稳定 ， 对 海南 岛 红 树 
VPD 同 SAL 相似 ， 对 红 树 林 GPP 的 


高 ; 此 外 ,由 了 


[树林 生态 系统 GPP 


I 6 一 1 


整体 上 呈现 出 旱季 


向 雨季 过 渡 期 的 


， 在 雨季 (6—10 H) ， 随 着 SST 的 逐渐 升 高 ，SST 成 为 主要 的 环境 
光合 速率 降低 可 能 与 SST 超过 红 树 最 适 SST 而 使 其 


酶 活性 降低 有 关 ; 此外， 


限制 因子 ， 


由 于 SST 


0 月 受 SST 的 限制 也 由 


pA 


Hz 
pA 


和 降水 分 布 的 
海南 岛 地 形 


民 制 很 小 ， 虽 然 温 度 较 高 会 使 得 
因而 VPD 对 红 树 林 GPP 的 限 
区 的 红 树 林 生 态 系统 GPP 表现 出 较 


西北 向 东南 
林 的 GPP 影响 很 小 。 
VPD 升 高 ， 但 伴随 雨 
制 低 。 从 整体 上 看 ,海南 岛 西部 、 北 部 


大 的 月 均 GPP 波动 ， 这 可 能 由 于 海南 岛 


部 为 


WA. 东部 为 迎风 坡 ， 受 东南 季风 影 
间 高 周边 低 ， 东 部 和 北部 多 为 丘陵 山地 ， 


海南 岛 北部 也 易 受 东南 季风 影响 。 相 较 之 下 , 海南 岛 南 部 受 五 指 山 阻 断 东 南 季风 ,气候 条 件 
较为 稳定 ， 包 括 相 对 较为 恒定 的 水 热 分 布 ， 因 而 月 均 GPP 值 相对 较为 稳定 。 


4 结论 


本 研究 利用 Sentinel-2 影像 和 环境 因 
年 海南 岛 红 树林 生态 系统 GPP 时 空 分 布 格 


局， 并 探讨 了 


形成 机 第 o 


= 


结果 表明 ， 


子 数 据 ， 基 于 MVP-LUE 模型 估算 了 2016—2020 


研究 期 间 


内 ， 海 南 岛 红 树林 生态 系统 年 GPP 呈现 出 增长 趋势 ， 在 空间 上 呈现 东部 高 于 西部 ， 北 部 高 

于 南部 的 分 布 格局 ， 在 形成 机 制 上 ， 海 南 岛 红 树林 生态 系统 月 均 GPP 时 空 分 布 格局 受到 多 

气象 因素 综合 影响 ， 年 内 季节 性 上 GPP 在 “ 旱 一 雨 ”过 滤 期 随 着 Tar 的 限制 降低 和 PAR. 与 

SST 的 促进 作用 ， 展 现 出 较 高 的 GPP; 早 季 中 期 Tir 低 ， 对 GPP 的 限制 严重 ， 但 随 着 纬度 

的 降低 ， 红 树林 生态 系统 GPP 受 Ta 限制 的 程度 被 减弱 ， 雨 季 则 由 于 较 高 的 云 覆 盖 ， 使 得 

PAR 成 为 GPP 的 限制 因素 。 
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